A Simple Model for the Deformation-Induced Relaxation 

of Glassy Polymers 

S. M. Fielding,! R. G. Larson,^ M. E. Cates^ 
^ Department of Physics, Durham University, Science Laboratories, South Road, Durham DHl 3LE, UK 
^ Department of Chemical Engineering, University of Michigan, Ann Arbor, MI 48109-2136, USA 
^ SUPA, School of Physics and Astronomy, The University of Edinburgh, 
JCMB Kings Buildings, Edinburgh, EH9 3JZ, UK 
(Dated: October 19, 2011) 

Glassy polymers show "strain hardening" : at constant extensional load, their flow first accelerates, 
then arrests. Recent experiments have found this to be accompanied by a striking and unexplained 
dip in the segmental relaxation time. Here we explain such behavior by combining a minimal model 
of flow-induced liquefaction of a glass, with a description of the stress carried by strained polymers, 
creating a non-factorable interplay between aging and strain-induced rejuvenation. Under constant 
load, Hquefaction of segmental motion permits strong flow that creates polymer-borne stress. This 
slows the deformation enough for the segmental modes to re- vitrify, causing strain hardening. 



PACS numbers: 64.70.pj,62.20.-x,83.80.Va 

Understanding the flow of polymeric materials is cen- 
tral to their manufacture and performance. After decades 
of progress, the flow properties of molten polymers are 
elegantly described by modern entanglement theories 
[T1I2]. In use, however, most polymeric materials are not 
molten, but rigid. This conversion is commonly achieved 
by cooling to below the glass transition temperature, Tg. 
In contrast to the molten case, satisfactory theories of 
polymer glass rheology remain elusive. 

Just below Tg, polymer glasses undergo slow plastic 
deformation if stress is applied [3l H] . Similar plasticity 
is shown also by molecular, metallic, and colloidal glasses 
[5H7]. Our understanding of flow in such non-polymeric 
glasses has improved greatly due to recent advances in 
microscopic [H [H] and mesoscopic [TUl - fTS] theory. Crucial 
to glass rheology is physical aging: a quiescent glass be- 
comes more sluggish with time, rejuvenating under flow. 
This is captured schematically in minimal 'fluidity' mod- 
els, with a time evolution equation for a single structural 
relaxation rate (the fluidity) FH] . In the so-called 
'simple aging' scenario, the structural relaxation time r 
(or inverse fluidity) of the system at rest increases lin- 
early with its age [H [TUl IS] • A slow steady flow cuts off 
this growth at the inverse flow rate. 

In polymeric glasses, new properties emerge from the 
interplay between polymeric and glassy degrees of free- 
dom. Particularly striking is the evolution of the seg- 
mental relaxation time T(t), controlling the rate of local 
rearrangements, when a load is applied. Recently, Lee 
et al [31 [15] showed that r(t) falls steadily during the 
early stages of elongational deformation, and then more 
sharply, reaching a small fraction lO""^'^ of its initial 
level before dramatically rising again, as the local strain 
rate started to drop on entering the 'strain hardening' 
regime. While elements of this scenario have been con- 
firmed in coarse-grained and molecular simulations [171 - 
[20] , no convincing theoretical picture has yet emerged. 



In [3], the results for r(t) were found inconsistent with 
the early theory of Eyring [21] and also with a more re- 
cent model [11] (see also [22H24] ) involving broadly sim- 
ilar precepts. Indeed the Eyring-like assumption of a 
purely stress-dependent fluidity, as introduced for poly- 
mers in |25j . is fundamentally at odds with aging in 
glasses, whose fluidity is time-dependent at constant 
stress [H [To] [16] . Previous work to incorporate aging and 
flow-rejuvenation into polymer glass theory has led to 
the Eindhoven Glassy Polymer model (EGP) [SH], where 
viscosity is controlled by a state parameter S that is 
age- and strain-dependent. However, in the EGP model 
treats aging and rejuvenation have factorable effects on 
S: strain-induced rejuvenation causes cumulative losses 
of structure (reductions in S) which multiplicatively re- 
duces all subsequent relaxation times. This is not what 
theories of simple glasses predict [IHl [H [H] . The EGP's 
precepts may thus be unsuited to the regime of strong 
fluidization, as addressed experimentally in [3] and in re- 
cent glass rheology theories [SHTUl [T^ [Tl] . 

Despite several recent efforts [TTl [22l [Ml [2Z] , creating 
a comprehensive theory of rheological aging in polymer 
glasses remains a formidable task. Here we show that a 
minimal model, combining just two key elements of any 
such theory (nonfactorable aging/rejuvenation, and the 
strain dependence of polymer-borne stresses), semiquan- 
titatively explains many of the results reported in [3]. 

Our model describes polymeric dumb-bells [51 sus- 
pended in a glassy 'solvent', whose microscopic relaxation 
time obeys a fluidity-type equation showing simple aging 
and flow-rejuvenation. Despite our nomenclature, we do 
not require any actual solvent to be present: the sepa- 
ration between polymer and 'solvent' instead divides the 
slow degrees of freedom of large sections of chain from 
the shorter-scale and faster relaxing modes that control 
local segmental dynamics. Our model thus follows lines 
developed in [221 US but crucially differs in its treat- 
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ment of aging and rejuvenation. For simplicity we treat 
the dumb-bells initially as purely elastic elements - as is 
valid in the molten state, where the elasticity is of en- 
tropic origin [2]. However, we later return to discuss the 
true nature of the polymer stress in polymeric glasses 
which is not solely entropic in character [231 HH] • 

Our model first defines a deviatoric stress tensor S = 
Gp{(tP - I) + Cia" - I) where ctP and cr" are dimen- 
sionless conformation tensors for polymer and 'solvent', 
G^'" associated elastic moduli (see below), and I the unit 
tensor. We then adopt the following equations for the 
conformation tensors and solvent relaxation time r: 

o-P + v.VcrP ^ crP.Vv+(Vv)^.(TP-a(crP-I)/T(l) 

&'+^r.\/(T' = c^^Vv+ (Vv)^.cr'' - (cr^ -I)/r (2) 

f + v.Vt = l-(T-ro)A (3) 

A(D) = ^V2Tr(D.D) (4) 

Here v is the fluid velocity and D = (Vv -I- (Vv)^)/2. 

Eq.Q is an upper-convected Maxwell model [2], de- 
scribing the dynamics of elastic dumb-bells; these carry 
an elastic stress G^a-P and have a structural relaxation 
time tP = r/a, proportional to, but much larger than, 
that of the 'solvent', r. In the simplest models of dense, 
molten, but unentangled polymers, a = N^'^ with N the 
polymerization index [2] , whereas in a lightly crosslinked 
elastomeric network [3] one expects a — 0. Consistent 
with its glassy nature, the solvent itself is viewed as a 
viscoelastic fluid. Bearing in mind that it represents 
shorter-scale polymeric degrees of freedom, we model this 
fluid using another upper-convected Maxwell model 
Because there are more local than chain-scale degrees of 
freedom, we expect G'^ > Gp . 

Finally, the solvent's structural relaxation time t obeys 
a fluidity- type equation (|3]), with the following two fea- 
tures. First, without flow, r increases linearly in time 
at a (dimensionless) solidification rate f(D = 0) which 
for simplicity we set to unity. This embodies the sim- 
ple aging scenario that emerges from mesoscopic models 
[lOj , whereby local configurations evolve into ever deeper 
traps. Second, with flow present, r would, in the absence 
of such aging, itself undergo deformation-induced relax- 
ation towards Tq which is a 'fully rejuvenated' value. This 
relaxation occurs at a rate A, proportional to a scalar 
measure of flow rate (with /i another dimensionless co- 
efficient mil]). In steady shear (A = /i-y), t then varies 
inversely with strain rate 7 in accord with microscopic 
theory [8]. For uniaxial elongation at strain rate e, Q 
reduces to A = /j,-\/3|e|. Note that in this simple fluidity 
model, the rejuvenation of r is essentially strain-induced 
[5] but, in contrast to the factorable model of [21], can 
be rapidly reversed by subsequent aging. 

Our model is completed by the standard equations of 
mass and force balance for an incompressible fluid of neg- 
ligible inertia: V.v = 0, and V.[I] -t- Irjn] = 0. (Here rj 
is a small additional Newtonian viscosity, included solely 



for numerical reasons [IH].) We have solved our model 
numerically for uniaxial extension flows within a lubrica- 
tion approximation appropriate to long cylindrical sam- 
ples. Our numerical solutions address two cases [2^ . One 
is an effectively infinite cylinder that remains of spatially 
uniform, but time dependent, cross section. The second 
addresses a finite cylinder perturbed to trigger an inher- 
ent 'necking' instability. The latter is commonplace in 
elongation of polymer glasses (and reported in [3] , in rel- 
atively mild form). We show next however that a semi- 
quantitative account of the r(i) response under elonga- 
tional load is already predicted by applying our simple 
model to the infinite uniform cylinder. 

In confronting the experimental data for r(t) we first 
set a negligibly small, appropriate for a crosslinked mate- 
rial [3]. The experimental protocol of [3| determines the 
applied tensile force F; the initial relaxation time {t^ in 
our model) before applying the load; and the time t„ at 
which unload later occurs. There remain four material 
parameters in the model: Gp,G^,to and fi. As detailed 
in three of these are strongly constrained by mea- 
surements that do not involve the dip in the r(i) curve. 
Indeed, Gp/F can be deduced from the asymptotic defor- 
mation in the strain-hardened regime just before unload; 
once GP is known (we find Gp = Q MPa) G"* and /z are in 
turn estimated from the step-change in r during initial 
loading, and from the separately measured slope [3] of the 
'effective flow curve' s(t). Hence the only unconstrained 
parameter in fitting the dip in T{t) is tq. 

We find a good semiquantitative account of the strain 
curve and r(i) data, up to but not beyond the point of 
unload, by choosing tq ~ 6s. (Unloading is addressed 
separately below.) Fig. [l] shows not only the local strain 
and the segmental relaxation time T{t), but also the ten- 
sile stresses TP''^ carried by polymer and solvent respec- 
tively. Key features of the experimental data, reproduced 
by our minimal model, include: (i) the initial drop in 
T on applying the load; (ii) its subsequent further de- 
cline to a state of strong fluidization, with a sharp min- 
imum Tmin near the point of maximum elongation rate; 
and (iii) its rapid rise from the minimum to a strain- 
hardened plateau prior to the point of unloading. Not 
only the initial tenfold drop in r on loading but also the 
subsequent further sharp dip is quantitatively accounted 
for. Figure [2] shows r function of the elongational 
stress, with breakdown of the Eyring-like expectation of 
a monotonic, single- valued plot. Figure [s] shows (on log- 
log) e against 1/t; this plot was found to collapse the 
experimental data in ^ and a similar, if lesser, effect is 
seen here. Considering the crudeness of our model (for 
example, the representation of polymers and solvent by 
a single mode each), this is remarkable agreement. 

If our model is correct, the physics of all these effects is 
remarkably simple. The (pre-aged) 'solvent' glass has a 
yield stress Ey (in our model this obeys Ey = G^g{y/3fi) 
with g{y) = 3?//(y — 2)(y + 1)) which is initially exceeded 
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FIG. 1: (Color online) Left Panels, solid curves: local strain 
e = expe — 1 |3], reduced relaxation time r(t)/tu, and tensile 
stresses T'''" = C'" {a^'f —a^'J) of the polymer (p) and solvent 
(s) during loading of an infinite uniform cylinder. Parameters 
Qs/QP ^ g 5^ ^ ^ ^2.5, t^/ro = 10^, TO = 6s; applied force 
/ initial area / — 2.7 G^. (The curve for , in red, initially 
lies below but crosses it during strain hardening.) The 
unload results for the basic model (6 = 1) is shown dashed; 
the solid curve after unload has 6 = 0.1. The horizontal axis 
is marked both in dimensionless model units (top) and real 
time (converted using tq), bottom. Right Panels: Compara- 
ble experimental data for local strain and reduced relaxation 
time. (From [3]. Reprinted by permission of AAAS.) 



by the applied load. After an initial step-down in r 
caused by step strain on loading, the material yields and 
progressively fluidizes further; accordingly its strain rate 
accelerates, giving positive feedback and a collapse in t. 
As deformation builds up, however, an ever growing share 
of the applied stress is instead carried by the stretching 
polymer chains. This causes the flow rate to drop, so that 
the solvent, whose stress now obeys < Sy, starts to 
solidify. This simple view of strain hardening also di- 
rectly explains the remarkable behavior of T(t). 

Models that factorize aging and rejuvenation effects 
are seriously challenged by the rapid recovery of r 
after the dip. (A multimode spectrum '30] is unlikely to 
help here.) With simple aging, such factorization predicts 
T ~ (t + iiu)/(e), so that if the segmental relaxation times 
falls from its pre-deformation value to a small value 
T — ftw = Tmin at the dip, a tenfold recovery to r ~ 
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FIG. 2: (Color online) Reduced relaxation time r/t^ 
against actual stress in the infinite uniform cylinder. Pa- 
rameter values C/G^ — 8.5, /i — 12.5, t^/ro = 
10*, To — 6s for a scaled applied force per initial area 
f/Gp = 0.3,0.6,0.9,1.2,1.5,1.8,2.1,2.4,2.5,2.6,2.65,2.7 (in- 
creasing left to right). The unload time obeys Tunioad = 
1500ro. The horizontal axis is marked both in dimension- 
less model units (top) and laboratory stress (bottom, as used 
in the inset; converted factor G^ = 6 MPa). These curves 
show qualitative agreement with the experimental data (from 
3^; reprinted by permission of AAAS) (inset). 



case of a finite cylinder subject to a necking instability. 
More details, and an additional figure, are provided in 
|29) . Although our model is not predictive of sample 
shapes (which depend on the details of the perturbation 
used to initiate the neck), plots of r(t), and sample radius 
p(t), at three different initial positions along the sample 
are in qualitative accord with the experiments of [Sj . The 
explanation given above for the temporal behavior of T{t) 
during elongation of an infinite uniform cylinder remains 
equally valid for a finite, necked one. 

In order to confirm that our model also behaves reason- 
ably in strain-controlled flows, we have calculated stress 
responses for startup of steady elongation and compres- 
sion. These show an overshoot (see additional figure in 
[29]) whose height varies as In(et^), similar to the be- 
haviour found in simple aging fluids [lOj, and in broad 
accord with the polymer glass literature. 

These successes are very encouraging. However, the 
model as formulated so far breaks down badly when the 
sample is unloaded. Here the experiments show a mod- 
est drop in r immediately on removing the load, followed 
by a gradual recovery towards the pre-deformation value. 
The dotted line in Fig. [l] shows the prediction based on 
Eqs. ([l]j4|; r drops, but then falls much further before re- 
covering. The reason for this behavior within our model 
is clear. In the strain-hardened regime, the polymers 
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FIG. 3: (Color online) For the same runs as in Figj2) during 
loading phase only, log-log plot of reduced strain rate against 
reduced relaxation rate. Inset: experimental data collapse 
with this plot, slope 0.92. (From [3]. Reprinted by permission 
of A A AS.) Partial collapse occurs here (dashed line is slope 
1): while re-entrant regions do not fully superpose, the slopes 
of the rising and decreasing curves remain comparable. 



carry a large elastic tensile stress, which exceeds Sy. 
Upon unloading, this acts backwards on the vitrified sol- 
vent, causing it to yield. The resulting t{C) resembles 
a re-run of the initial loading experiment. Another dis- 
crepancy is that the value of G'' ~ 6 MPa needed to fit 
the loading data is approximately ten times larger than 
the rubbery modulus of the same material above its glass 
transition (see, e.g., [21]). This confirms that the strain- 
hardened modulus of polymer glasses does not primarily 
stem from single-chain entropic elasticity [211 [IH] . 

We now identify a physical mechanism that could ac- 
count for both discrepancies. We invoke the well es- 
tablished phenomenon of large but viscous stresses that 
arise when chains are strained rapidly relative to their 
own relaxation time (er^ ^ 1). Under such conditions, 
relatively small sections of the polymer quickly stretch 
close to full extension locally, forming a nearly one- 
dimensional multiply folded ('kinked') filament [5^ I33j . 
Further stretching occurs by migration and annihilation 
of neighboring kinks of opposite sign. During this pro- 
cess, a large fraction of the stress carried by the polymers 
is not entropic-elastic, but instead caused by viscous drag 
against extended subsections of chain. Upon unloading, 
a large fraction of this inelastic polymer stress disappears 
on a very rapid timescale |32j . This mechanism is closely 
related to the phenomenon of 'chain conformation hys- 
teresis' which arises for fully stretched chains, causing 
a sudden loss of polymer stress on unloading with only 
modest relaxation of polymer conformations [5^ . 

Rather than attempt a full treatment of this rather 
complex effect (which would certainly require a multi- 



mode polymer description |32|). we retain our equations 
but suppose phenomenologically that the effective poly- 
mer modulus drops by a certain factor, QG^ , dur- 
ing unloading of the sample. The solid line in Fig. [T] 
shows the result for B = 0.1. This choice of B is consis- 
tent with the fitted being ten times larger than the 
value expected from entropic elasticity alone. The poly- 
mer stress acting backwards on the solvent is now safely 
below the solvent yield stress; the result is a modest drop 
and then slow increase in T(t), as seen experimentally. 

Overall, the success of our simplified model suggests 
that the striking time dependence of the segmental mobil- 
ity under elongation, reported in [3], should be a robustly 
universal feature of near-Tg polymer glasses. However the 
quantitative details strongly depend on dimensionless pa- 
rameters such as G'*/G^,/x and B. We cannot link these 
directly to microscopic physics, but such parameters can 
influenced by increasing polymer stiffness, adding small 
molecules, or introducing short side chains. (All of these 
should increase G^'/G'', by raising the ratio of solvent-like 
to polymeric degrees of freedom.) Our model may thus 
suggest design strategies for manipulating the evolution 
of t(<), tailoring the mechanical responses of polymer 
glasses to suit particular design needs. 

In conclusion, we have presented a simple approach to 
the modeling of polymer glasses that builds on recent 
models of rheological aging and rejuvenation in simple 
glassy fluids. Without attempting to capture every fea- 
ture of the experiments of [3] (for instance, we do not 
address the non-exponential form of local relaxations), 
the minimal combination of a simple-aging fluid with 
a strain-dependent polymer stress can explain much of 
what happens when a polymer glass is subjected to elon- 
gational load. The unloading behavior is less easily ex- 
plained, but consistent with a plausible modification of 
the same model, which crudely allows for the presence of 
non-elastic polymer stresses when er^ is large [35] . 

Our work suggests that an accurate representation of 
aging and rejuvenation physics will form a key part of 
any more comprehensive theory of polymer glass rheol- 
ogy. It encourages the view that a more comprehensive 
account of polymer glasses might be achieved by judi- 
ciously combining existing types of nonlinear rheologi- 
cal theory (describing non-glassy polymers and simple 
glasses respectively). Quantitative progress along these 
lines might enable rapid advances towards the design of 
superior polymer glass materials. 
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SUPPLEMENTARY MATERIAL 

Here we detail the model equations and their trans- 
formation via a lubrication approximation into a form 
suitable for numerical study (by standard methods that 
we mention only briefly) . This is followed by presentation 
of additional results on necking and on startup flows, and 
a discussion of the parameter estimates made in relation 
to the experiments of Lee et al [3] . 

Equations of Motion and Lubrication Approximation 

We adopt Equations (1-3) of the main text plus the 
continuity equation for an incompressible fluid (V.v = 0) 
and the force balance equation which reads 

V.[S + 27?D]=0 (1) 

with the deviatoric stress written as 

1:^9GP{(tP + G' (cr" - I) (2) 

Here the phenomenological factor 9 is unity during load- 
ing but can drop upon unload as described in the main 
text. The parameter 77 is a small additional viscosity in- 
troduced for purely numerical purposes (to avoid having 
to deal with inertia), which is always negligible in prac- 
tice. 

To constuct a lubrication approximation (long slender 
samples) we introduce a coordinate z along the elonga- 
tion direction; we denote by x a radial coordinate and 
A{z,t) the cross sectional area. The local fluid velocity 
along z is denoted v{z, t) so the local elongation rate is 
e — dzV. The continuity equation then reads 

A = -vd,A - Ae (3) 

whereas force balance demands 

dzF = (4) 

where the local tensile force is denoted as 

F{t) = [^(S,, - S,, + 377e)] (5) 

(Force balance require this to be independent of z.) In 
this expression for F we have used the lubrication approx- 
imation that velocity gradients depend on z only (not x) 
from which (Vv)^^ = —2{'Vv)^x follows by incompress- 
ibility [55H37j . The same approximation in the Maxwell 
equations for polymer and solvent (Eqs. 1 and 2 in the 
main text) gives 
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with Tp = r/a the polymer relaxation time. The equa- 
tion of motion for t (combining Eqs. 3 and 4 of the main 
text) becomes in the lubrication approximation 

f = -wa^r + 1 - (t - ro)Af\/3|e| (10) 

Coordinate Transformation 

We now specialize to the case where the applied tensile 
force F(t) is a piecewise constant function of time. We 
make a coordinate transformation that removes the expo- 
nential increase in sample length (and matching shrink- 
age of area) corresponding to uniform afline deformation 
given by the spatially averaged strain e{t). This greatly 
helps the numerical analysis by obviating the need for 
an adaptive numerical mesh to cope with large sample 
deformations at late times. The variables z, w, e, A, F are 
transformed into a new set m, w, C, a, / as follows: 



z 


= uexp[e(i)] 


(11) 


V 


= \e{t)u + w{u, t)] exp[£(t)] 


(12) 
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= j(t) + C(u,i) 


(13) 
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= Aoa(u, i) exp[— e(t)] 


(14) 


F 




(15) 



where / is (piecewise) constant in time. The remain- 
ing variables (c^^ cr^^ r) are not transformed, but 
expressed as functions of (u, t) rather than (z, t) as pre- 
viously. Without loss of generality we take the initial 
length of the sample to be unity; in the absence of neck- 
ing, the sample domain is < m < 1 for all subsequent 
times. The initial cross sectional area ^0 defines the un- 
deformed starting shape of the sample, with a{u, 0) = 1 
for all u. 

The transformed equations of motion are 
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To solve these equations numerically, zz^'^'xx zzt'^ 

are updated each time step using Equations [T7p2| Up- 
dated values are fed into Equation [THj whose integral 
over the entire sample determines e; this value is then 
substituted back to determine C,(u^£) everywhere. 

If one assumes the sample does not neck (so that the 
deformation remains affine) these equations reduce to or- 
dinary differential equations: quantities depend on t but 
not w, and by construction, C = 0. Numerical solution 
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FIG. 1: (Color online) Left Panels: the reduced sample ra- 
dius p/po (with po the initial value), the local strain measure 
e — expe — 1 [3^, and the relaxation time T{t) during sample 
loading for three different positions in a finite cylinder subject 
to necking instability. Parameters as in Fig.f of main text, 
now for a cylinder that initially occupies the space z = to 
z = 1, spatially resolved using 400 numerical grid points. The 
three positions chosen are initially at z = 0.5,0.9,0.95. Right 
Panels: Similar experimental data; note the discrepancy in 
horizontal scales (no such discrepancy arose for Figs. 1-3 of the 
main text). (From [3]. Reprinted by permission of AAAS.) 
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FIG. 2: Upper Panels: Dimensionless tensile stress {T^ + 
T^)/G^ versus Hencky strain it in startup of steady elonga- 
tion at strain rate e. Left I^/tq — 10* with (bottom to top) 
ero = 10-^10-^10-'',10-^10-^ Right, ero = 10"^ with 
(bottom to top) t„/ro = 10^,10^,10^. Lower panels: The 
same, for startup of compressional flow (e < 0). The remain- 
ing model parameters are as in Fig.l of main text (except that 
the (unimportant) numerical viscosity rj is set equal to zero). 



proceeds by an explicit Euler algorithm [55] , with careful 
attention paid to timestep convergence. In our numer- 
ics we chose a — 10^^^, but the results are robust to 
variations at small a. We choose t] = 0.05, which, as pre- 
viously indicated is effectively zero and the same applies. 



Necking 

In general, however, one expects non-afhne deforma- 
tion of the sample due to necking. The partial differential 
equations (16-221 must then be solved in full. For this 



we split the operator, using explicit Euler for the local 
terms and first order upwinding for the convective ones 
[38j . In addition, it is now necessary to specify bound- 
ary conditions, and also to introduce some perturbation 
to the affine state so as to trigger the necking instabil- 
ity. In practice, necking occurs not in a random position 
but near the middle of the sample - presumably because 
clamping or other boundary conditions at the extrem- 
ities break the translational symmetry. Retaining pe- 
riodic boundary conditions for numerical simplicity, we 



add an early time perturbation of the form 



= [exp(-M^exp(2£)/Ao)+exp(-(l-M)^exp(2e)/Ao)] 



(23) 

whose effect is to create a sample comprising a long cylin- 
drical section that flares out at both ends on a length 
scale ~ \fAQ comparable to the sample width. (This is 
a plausible representation of the early-time shape of a 
sample that is not yet necked but whose ends, for what- 
ever reason, have an area different from the affine value.) 
This perturbation is applied as an additional term in Eq. 
[rf| but only at early times; it is switched off once the 
mean strains e exceeds a chosen threshold, e* . 

Fig. [T] shows results for the sample radius and seg- 
mental relaxation time at three different positions in a 
necking cylinder, computed within the lubrication the- 
ory described above. 



Strain-Controlled Flows 

Figure [2] shows our predictions (for a non- necking infi- 
nite cylinder) for the stress overshoot in polymer glasses 
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in startup flow. We are not aware of startup data on 
the corresponding experimental system as studied by [3] . 
Nonetheless these curves qualitatively resemble literature 
data on startup in polymer glasses (see, e.g., Fig. 1 of 
PB]). Our model predicts a logarithmic dependence of 
overshoot amplitude on both strain rate and sample age 
with roughly a two-fold overshoot increase for a ten-fold 
increase in age or strain rate. (This ratio may depend on 
model details.) The overshoot is controlled mainly by the 
underlying fluidity model rather than polymeric effects. 



Parameter Estimation 

As mentioned in the main text we set a ~ 10^^^, which 
is effectively zero, and set the dimensionless aging param- 
eter dr/dt (in a system at rest) to be unity for simplicity. 
We have already chosen our length unit as the initial sam- 
ple length. The parameters remaining in the model are 
the material parameters G^, G*, tq, /j.; the tensile force F; 
the waiting time t^; the unloading factor 6; the pertur- 
bation threshold e* introduced above; the initial cross 
sectional area of the sample Aq] and the additional vis- 
cosity TJ. 

As stated previously, the viscosity ij is a numerical 
device that allows force balance to be maintained from 
one time step to the next by slight adjustments of the 
strain rate. (An alternative would be to include fluid in- 
ertia; another would be to have an extra iteration loop to 
maintain precise force balance at every timestep.) Any 
small value for rj will achieve this without corrupting the 
physics of the model; we choose ij = O.OSG^tq. Since 
C > and r > tq, this always remains negligible com- 
pared to the solvent viscosity. 

The initial cross section Aq can in principle be matched 
to the samples reported by Lee et al [3], but in a homo- 
geneous deformation none of the quantities we report de- 
pend on sample shape. Hence Aq is relevant only when 
necking is accounted for; when a definite value is needed 

0.1. (Recall that 



(as in Eq 23 above) we choose 
the unit of length is the sample length) . Similar remarks 
apply to £*; for definiteness we choose e* — 0.04. 

The applied tensile force F = AqJ can be estimated 
by hypothesizing that, in the plateau that is reached af- 
ter strain hardening, this force is carried predominantly 
by the polymers. Thus F = A{t)GP{aP, - <jP^) ~ 
A(t)GPexp(2e); it follows that / GPexp(e). Lee et 
al [3] report deformations in terms of a local strain vari- 
able {L{t) — Lq)/Lq which equates to exp[e(i)] — 1 in 
our notation. Mildly extrapolating the data of their 
Fig.l to estimate the plateau height, we conclude that 
exp[£(t)] ~ 2.7; it follows that / = 2.7GP. 

Note that / = F/Aq is in fact the 'engineering stress' 
of which is 16 MPa in the loading experiment re- 
ported in their Fig.l. The resulting G^ of ~ 6 MPa is 
about ten times larger than the rubbery modulus of the 



same material [3T]. This is consistent with the view of 
[24] that the modulus responsible for strain hardening, 
though attributable to polymers, is not simply entropic 
in origin as in rubbers (see main text). 

We turn next to the parameter ^ = X/{^/3\e\) which 
controls the strain-induced relaxation of r. Within our 
model, this is fixed by observing that in steady state at 
low elongation rates (|£|to ^ 1) the solution to Eq|lO| 
approaches 



(24) 



Figure 3 of Lee et al [3] plots £r(0) against t{0)/tss (with 
r(0) ^ To the pre-deformation value of the relaxation 
time); the slope on this log-log plot is close to unity. 
Estimating the near-constant ratio s/t~^ from the left 
edge of the experimental plot gives e — 0.045/t5s from 
which we deduce that ~ 12.5. 

Armed with this value, the solvent modulus G* is now 
estimated from the virtually instantaneous drop from 
the undeformed value of r = r(0) to a smaller value 
T+ ~ 0.1t(0), reported immediately after imposition of 
the load. The resulting strain £"*" is not directly reported 
in |3] but by integrating our equations of motion is found 
to obey 

/exp[£+] = [GP + G'')(exp[2£+] - exp[-£+]) (25) 

so that for small strain and large G^/Gp we have 
G'* = //(3£+). Integrating the r equation through 
the same step strain, we find the post-load value t'^ — 
T(0)exp[-V3Ai£+] = r(0)exp[-//;i/V^G^]. Figure 1 in 
[3] shows t+/t(0) ~ 0.1 from which it follows that 
G7GP = 8.5 (using the previous result that / = 2.7GP). 

The resulting G" ^ 50 MPa is 30 — 50 times smaller 
than the Young's modulus of fully vitrified PMMA [55] . 
which suggests that in the temperature range of interest 
(near but below the glass transition) only a small pro- 
portion of the monomeric degrees of freedom are directly 
involved in the process of solvent arrest. 

The next parameters to be estimated are tw , the effec- 
tive age of the sample on first loading, and tq which is the 
relaxation time of a fully devitrified sample. Within our 
model, the initial relaxation time r(0) is simply equal to 
t^. The remaining parameter tq is then needed to con- 
nect the time axis in our numerical calculations to that 
in the experiments of [3]. In essence, we use this as an 
implicit fit parameter to get the best agreement between 
our r(t) curves and those reported experimentally. (In 
numerical practice, we use tq as the unknown but fixed 
unit of time in which case t^/TQ, the numerical sample 
age, is the explicit fit parameter.) To extract a value of 
Tq from our fits we observe that, when the loaded sample 
in Figure 1 of [3] achieves its minimum relaxation time 
(point (e) on panel C of that figure), the decay curve 
from which this time is measured (panel B) is close to 
mono-exponential with a relaxation time of 28s. The re- 
laxation time t(0) for the undeformed material (for which 
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the decay is far from exponential but fitted carefully by 
the authors of [3]) is (from their Figure 1 A) 10^ "^ times 
larger, i.e., 6 x 10^ s. Our best fit waiting time obeys 
tw/'To = 10^; since t(0) — t^, we then have tq = 6 s. 

As noted in the main text, the elongational dynamic 
yield stress of the solvent fluid can be found by calcu- 
lating the steady state stresses at fixed e and studying 
the limiting tensile stress T as e — ^ 0. The result is 
Sy = \/3G*//i ~ 7 MPa. This is somewhat at odds 



with Figure 2 of [3], which suggests an elastic response 
up to stresses of about twice this value. However, the 
true yield stress could lie below the value extracted from 
the experimental measurements if the latter do not de- 
tect extremely slow creeping motion that may arise for 
stresses close to but above Ey. 

The above considerations fix all the parameters of the 
model except for the unloading factor 9 which was ad- 
dressed already in the main text. 



